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We have numerically investigated whether or not a mean-field theory of spin textures generate 
fictitious fiux in the doped two dimensional t — J-model. First we consider the properties of uniform 
systems and then we extend the investigation to include models of striped phases where a fictitious 
flux is generated in the domain wall providing a possible source for lowering the kinetic energy of 
the holes. We have compared the energetics of uniform systems with stripes directed along the 
(10)- and (ll)-directions of the lattice, finding that phase-separation generically turns out to be 
energetically favorable. In addition to the numerical calculations, we present topological arguments 
relating flux and staggered flux to geometric properties of the spin textureJj The calculation is based 
on a projection of the electron operators of the t — J model into a spin texture with spinless fermions. 



PACS number(s); 75.10.Jm 
I. INTRODUCTION 



It is well known that topological spiaptextures are im- 
portant in quantum Hall ferromagnets.0i3 It is also known 
that an electromagnetic flux through a system of tight 
binding electrons on a lattice can lower the electronic ki- 
netic energy .□ Topological arguments also suggest that a 
possible source of this flux could be the formation of a 
spin-texture, which would suggest a relationship between 
doping and magnetic flux.Q However, for the two dimen- 
sional jantiferromagnetic Heisenberg model it has been ar- 
guedH'Ll that, taking the continuum limit and looking at 
long wavelength fluctuations about the Neel state, there 
is no topological term in the effective action. Although 
this argument is correct, it does not, however, answer the 
question whether or not spin-textures can be important 
on a length scale comparable with the lattice. Further- 
more it does not address the issue of second-neighbor 
hopping. These ideas will be explored in the., present 
paper which extends ideas presented previously.EJ In par- 
ticular, we generalize the topological arguments given in 
Ref. [ y and present a numerical comparison of the en- 
ergy of flux generating spin textures and flux free spin 
configurations for uniform systems. Besides uniform sys- 
tems, we also consider if spin- and charge-stripes arise 
naturally as a topological fictitious flux generating spin 
texture. 

The paper is organized as follows: In Section we in- 
troduce the t — J-model and provide some background 
material on its properties. We proceed in Section III 
by deriving the effective model which we will work with. 
This model turns out to inc lude topological fluxes which 
are discussed in Section IV and in Section ^ we review 
the effect of such a flux on a system of free electrons. Sec- 
tion 



VI contains a numerical mean-field investigation of 



the energetics of the system, comparing flux generating 
states and more regular spi rallin g states where there is no 
flux generated. In Section VII we extend the discussion 



of uniform systems to include stripes forming antiphase 
domain walls between Neel ordered regions. We describe 
of stripe model and present data from numerical calcula- 
tions comparing different stripes. Finally, we summarize 



our results in Section VIII 



II. THE T J MODEL 

In order to explore these ideas we employ the t — J 
model which most simply captures the competition be- 
tween Heisenberg exchange and kinetic energy: 
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The summation is restricted to nearest-neighbor pairs 
on the square lattice and the spin operator is given by 
Sr = ^cJ„iT"^Cr/3, whcrc a = (cTx, CTj,, (T^) is the vector 
of Pauli matrices. All states containing doubly occupied 
sites have been excluded from the Hilbert space, leav- 
ing three states per site: | ), | j), and | [). A natural 
generalization, deferred to later in this paper, is to add 
Coulomb-repulsion between particles occupying nearest- 
neighbor sites. 

Striped phases have been found experimentally in high- 
Tc materials to which the present model has been ap- 
plied. There is an ongoing debate regarding the existence 
of stripes in the t — J- and Hubbard models. Stripes 
were first faund in Hartree-Fock solutions of the Hub- 
bard modelQ, but the stripes found in these calculations 
had one hole per unit length of the stripe, in contrast 
to the results from experiments where half a hole per 
unit stripe length is found. From DMRG calculations 
on finite .avstems (of the order 20 x 10). White and 
ScalapinoElilJ find stripes in a wide range of dopings. For 
instance, using J = 0.35t they find stripes for dopings in 
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the interval < x < 0.3. For x < 0.125 the stripes 
have half a hole per unit cell of stripe, in agreement with 
experiments, and the distance between two consecutive 
stripes is d = {2x)~^. For higher dopings they find that 
there is one hole per unit cell of the stripe and that 
the inter-distance between the stripes is d = x^^ . On 
the other hand, using quanj|iim Monte-Carlo calculations 
Hellberg and ManousakisOEj find that uniform or phase- 
separated states are energetically favorable. In this case, 
the formation of stripes would also require the existence 
of a long-range Couloynh i.nteraction preventing an ordi- 
nary phase-separationoEj. 

Incommensurate states were discussed in connection to 
the t — J model before the notion of stripas .was intro- 
duced. It was found by Shraiman and SiggiaE^EZI using a 
continuum limit of the model, that the antiferromagnetic 
order of the undoped t — J model is unstable against the 
formation of a spiral state for small dopings. JJsing var- 
ious mean-field approaches other authorsH^Lj came to 
similar conclusions. We will return to the spiral instabil- 
ity in Section VI using the effective model to be described 
in the next section. It will be shown how a small twist in 
the spin-order leads to a reduction of the kinetic energy 
of the order tx while the loss in exchange energy is of the 
order Jx^] showing that for small enough dopings there 
is energy to be gained by twisting the spin-order. 

In addition to spin textures, Affleck and MarstonEil 
discussed the possibility of flux textures. They replaced 
the two spin components of the electrons by N different 
flavors, extending the SU(2) spin symmetry to SU(A^). 
Taking the limit N oo, they obtained an essentially 
exact mean-field model to which they numerically looked 
for solutions. In particular they found a phase, called 
the flux phase, where the sum of the phases of the link- 
operators Xrr' = cjcr' around a plaquette equals ±7r. 
This is interpreted as half a flux quantum penetrating 
each plaquette. These phases do not come from a real 
electromagnetic field and are therefore referred to as fic- 
titious. |-| 

Work by Hasegawa et alB showed that the energy of 
non-interacting spinless fermions has a minimum when 
a uniform flux, corresponding to one flux quantum per 
particle, threads the system. States having ^ = n (in 
units of the flux quantum) are referred to as commen- 
surate flux states. This shows that a fictitious flux can 
lower the kinetic energy of the particles. The commensu- 
rate flux states have also been considered in connection 
with the t — J model. 

Another possibility for a flux state is to have a stag- 
gered flux through the system. In the case of half a flux 
quantum per plaquette there is no difference between 
uniform- and staggered fluxes, so the Affleck-Marston 
state can be thought to belong to this category as well. 
Insaired by the work of Shraiman and Siggia, Kane et 
alx2 suggested a double spiral state showing a staggered 
chiral spin opder, and hence also, according to a result 
of Wen et alB, a staggered fictitious flux. In a staggered 
flux state, the time-reversal symmetry can be broken lo- 



cally but not globally, as the system is invariant under 
a time-reversal operation followed by a lattice transla- 
tion, just hke a Neel state. Staggered jfijoi-phases have 
also been investigated by other group^ScJ. The effec- 
tive model used by Barford and Lrc3 coincides with the 
model derived in the next section. 

Our paper expands on results found by previous au- 
thors; we find that certain spin textures and charged 
stripes in particular are coupled by the creation of a fic- 
titious TT flux which we show is a natural consequence of 
a stripe with broken rotational symmetry. 



III. DERIVING AN EFFECTIVE MODEL 

In order to make progress, we make certain simplifi.- 
cations of the t — J model. First, following Schulzta, 
we introduce a local quantization axis llr at site r. 
In terms of spherical coordinates we write f2r = 
(sin 6'r cos (/)r, sin^r sin 0r, cos 6'r). This local SU(2) trans- 
formation on Cr, denoted by U-^ , must fulfill the equation 



(2) 



As can be seen from the above equation, specifying Sir 
determines J7jj only up to a rotation about the new local 

z-axis. For example, we may choose our SU(2) transfor- 
mation according to 
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where uj^ — Cz x D,,.) ~ sin 0r, cos 0r, 0). 

Expressing the t — J Hamiltonian in terms of this new 
spin-coordinate system we find 



^ = E [~t{cLM:icr^p+h.c. 
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with 



cos66ij + (1 — cos6)u!iU!j 
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We note that R-^ is the SO (3) rotation-operator induced 
by the SU(2)-transformation U-^. 

Thinking of the t — J model as the large-C/ limit of the 
Hubbard model, we know that there is a gap between the 
Hubbard bands scaling as U , corresponding to-the energy 
cost for a double occupancy. Following Schulaij, who ne- 
glected holes in the lower Hubbard band, we will throw 
away the upper Hubbard band because of this large gap 
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when we consider hole-doping, corresponding to the re- 
moval of states containing double occupancies from the 
Hilbcrt space of the Hubbard model. Since the quanti- 
zation axis at a site is locally determined by fir, we can 
arbitrarily assume that the upper Hubbard band is as- 
sociated with spin down relative fi. Hence, our effective 
model is obtained by keeping only the terms in Eq. (^) 
associated with spin-up particles. The spin of an electron 
at site r will now be determined by the field fir- As the 
simplest approximation, we will consider the fi-field as a 
classical field, neglecting spin-fluctuations in the system. 
Keeping only terms containing particles aligned with the 
positive local z-axis we obtain an effective Hamiltonian 



E 
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with t""' = tMll, K"''' = y{hr ■ f^r' - 1), and Cr = CrT. 
This Hamiltonian describes a system of spinless fermions 
moving in a lattice with position dependent hopping am- 
plitudes and interaction strcngths.EZI 

The effects of the spin texture on the charge-motion 
is therefore to generate an effective hopping amplitude, 
t — > t cos I and the appearance of a fictitious magnetic 
field. We also note that there is a Coulomb like nearest- 
neighbor interaction of strength ■|j(cos6' — f). When 
cos 6* < 1 this leads to an effective attraction between 
particles, which hints that the system may favor a phase- 
separation when being doped. 

We have already mentioned that there is a degree of 
freedom in not being fixed by fi. Since the effect of 
a spin-rotation about the local z-axis on the "up" spin 
only introduces a phase factor, it will be indistinguish- 
able from a local electromagnetic gauge transformation 
in our approximation. Hence, the set of physically in- 
equivalent choices of J7jj are determined by Q, i.e. they 
belong to SU(2)/U(1);^ S'^. 




FIG. 1. The mapping of a real-space plaquette into spin 
space. The solid angle spanned by the mapping is given by 

The physical degrees of freedom of the hopping are con- 
tained within the size of the hopping amplitude and the 
gauge invariant parts of the complex phases of the hop- 
ping elements. In case of nearest-neighbor hopping only, 
the smallest closed loop that can be formed is around a 



plaquette in the lattice, see Fig. and the flux enclosed 
by such a counter-clockwise path 0— >1^2^3^0is 
given by 
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One may show that this flux $0123 is equal to half the 
solid angle enclosed by the shortest path on the spin 
sphere connecting the points {riijf^g. Thus the flux is 
equal to 27r(5, where Q is the topological charge repre- 
sented by the plaquette. 

In the following discussion we will refer to the flux as 
a fictitious flux, in contrast to a "real" electromagnetic 
flux that would come from an applied magnetic field. The 
reason for this distinction is that the fictitious spin gen- 
erated flux is only seen by the spinless fermions in the 
system and is unrelated to a physical electromagnetic 
field. Furthermore, the fictitious flux does not couple to 
the charge of the fermions, but rather to the z-component 
of the spin measured in the local spin-coordinate system. 
Since all particles in our system are polarized along the 
positive z-axis, they will appear as having the same fic- 
titious charge. However, the flux can still drive currents 
through the system and in principle it is possible for the 
fictitious fiux to cancel the effect of an external electro- 
magnetic flux on the particles in the system. 

A physical magnetic field giving one flux quantum 
per plaquette is enormous. If we assume that the lat- 
tice constant is a ~ 5 A, the resulting magnetic field is 
B = h/{ea^) ~ 10^ T. This energy is much larger than 
the typical electronic energies per site, being of the order 
a few eV. 



IV. PROPERTIES OF FICTITIOUS FLUXES 

In Ref. [ P we investigated the relations between 
the magnitude of the hopping t^^ and its complex 
phases for both ferromagnetic and antiferromagnetic 
spin-configurations. In the present section we shall 
briefly review and generalize these arguments. 

We again consider a square lattice with a set of spins 
placed at each lattice site r. The interior angles on the 
surface of the sphere are described by angles as can be 
seen in Fig. |^. The fictitious flux through the plaquette 
is equal to half the solid angle covered by the plaque- 
tte in spin space, which by spherical geometry is given 
by the sum of the interior angles in excess of 27r. It is 
obvious from Fig. |I| that if Orr' , the angle subtended by 
the arcs on the sphere is small, the area of the spherical 
parallelogram, and hence the fictitious flux, will be small 
as well. The following expressions give the size of the 
hopping and the fictitious flux through the plaquette: 
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For small values of opening angles on the spin sphere 
6'rr' — and ~ a2+i the flux is approximately given by 
^0123 — — cos "^2"^ ' showing that the flux is bounded 
by|$|<f. 

If we instead turn our attention to antiferromagnetic 
configurations described by 6* > ^ , the situation changes 
drastically. The path taken in spin space when going 
around a plaquette in an antiferromagnetic configura- 
tion is shown in Fig. ^j. If we denote the antiferromag- 
netic (staggered) spin at site r by f2r, the path taken 

is r^o — * f^i — ^ ^^2 — ^ — > ^^Oj where Q! = —fl 
denotes the antipodal point on the sphere. 




FIG. 2. This figure shows the path taken on the spin 
sphere when going around a plaquette in an antiferromag- 
netic spin-configuration. 

Redefining 9 to denote the opening angle between two 
neighboring staggered spins, we find the following rela- 
tionships: 

Ir'-'-'l =tsin- 

*01'23' =^-^5I(-l)'"^ (9) 

i 

We note that the size of the fictitious flux is now com- 
pletely decoupled from the opening angle 9, as long as 
9 < ^. Furthermore, in the antiferromagnetic case the 
fictitious flux is staggered since the path on a neighbor- 
ing plaquette will be traversed with the sublattices ex- 
changed. 

If we allow for the possibility of next-nearest-neighbor 
(diagonal) hopping, we have to take into account addi- 
tional gauge invariant fluxes. There are four of these for 
each square plaquette, defined by the removal of one of 
the four vertices of the plaquette. In Ref. [ Q we show 
that there is a topological constraint relating these four 
fiuxes. This constraint takes the form 

'I'01'2 + $023' - $01'3' - $1'23' = 27m, (10) 

where n = for a ferromagnet and n = ±1 for an anti- 
ferromagnet. In case of a ferromagnet a prime does not 
denote the antipodal point, but rather the point itself. 



This relation is easy to verify by looking at Fig. y noting 
that the sphere is exactly covered by the four regions in 
Eq. (|lO|). The topological constraint, Eq. (|lO|), does not 
rely on the assumption of a spherical parallelogram and 
does also hold in the presence of an external electromag- 
netic flux. Counting degrees of freedom, we know that 
there are two degrees of freedom per plaquette (or site) in 
choosing the spin-conflguration and in addition we have 
one parameter coming from an external flux. All in all, 
there are three free parameters per plaquette and hence 
we expect that the four fluxes through the sub-triangles 
are related by a single constraint, given above. 

V. SPINLESS FREE FERMIONS WITH FLUX 

Before turning to a more through investigation of the 
physics of the effective Hamiltonian, Eq. (^, we will dis- 
cuss the effects of a flux through a system of free spin- 
less fermions. Hasegawa et alu investigated a system of 
free electrons on a square lattice with a uniform magnetic 
flux. They found that the energy is minimized when there 
is exactly one flux quantum per particle, i.e. the optimal 
flux per plaquette is related to the doping according to 
<I> = (1 — x)^o, where $o is the flux quantum. To il- 
lustrate this effect we have diagonalized the Hamiltonian 
for such a system of free fermions in a uniform flux $ for 
different values of finding the single-particle energies 
ei(<I>). The total energy of the system is then found by 
summing up the single-particle energies according to 

ei(<S>)<eF{7i,<S>) 

where epiji, $) is the Fermi energy corresponding to fill- 
ing n and flux $. 
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FIG. 3. This figure shows energy versus filling for free 
electrons using a set of different uniform fluxes, "1> = qn/12, 
where g = 0, 1, . . . , 12. In the center of the figure, i.e. at 
n = 1/2, the optimal energy is given by "l" = tt. Moving to 
the right, the next minima corresponds to $ = ll7r/12 and 
so on, finally finding an optimal flux $ = at n = 1. 
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In Fig. |[ we plot E{^, n) — E{0, n), i.e. the energy per 
site for different fluxes compared to the flux free case. 
The figure clearly shows how the optimal flux changes 
with doping. We also note that the system is particle- 
hole symmetric and hence E{^, n) — E{^, 1 — n). 

However, the flux that is generated by the antiferro- 
magnetic skyrmions is staggered in which case the Hamil- 
tonian can be exactly diagonalized. The spectrum of this 
system, assuming 0^ — Oy = 0, is 



Ht-t' 



ek($)=±2tsin( - 



X Y cos^ kx + 2 cos — cos kx cos ky + cos^ ky. (12) 

In Fig. ^ we show a plot of the energy per site for different 
staggered fluxes. The figure shows that the optimal flux 
is either or ±tt depending on doping, the only choices 
which are consistent with time-reversal invariance. 
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FIG. 4. This figure shows energy versus filling for free 
electrons compared to the fiux-free case for a set of different 
staggered fluxes, $ — qn. Going upwards in the center of the 
figure q — 1, 1/2, 1/3, 1/6. Note that for a certain doping, the 
minimum energy is obtained either for $ = or $ = ±7r. The 
crossing points occur at n ~ 1/2 ± 0.165. 

From this analysis we conclude that it is reasonable 
when searching for minimum energy spin textures to con- 
sider configurations supporting flux (coplanar configu- 
rations) and ±7r. In the next section we will construct a 
mean-field theory based on these observations. 



A. Second- neighbor hopping 



As we have seen in Section III, an antiferromag- 
netic spin-configuration strongly suppresses the effective 
nearest-neighbor hopping on the square lattice. However, 
a next-nearest-neighbor hopping is compatible with an- 
tiferromagnetic order. For the purpose of illustration let 
us consider the following purely kinetic Hamiltonian de- 
scribing spinless fermions 
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where ((rr")) denotes next-nearest-neighbor pairs. We 
remark that the sign of the nearest-neighbor hopping t 
is irrelevant as it can be changed by the transformation 
Cr I — > (— l)'"Cr. This transformation leaves the sign of t' 
unchanged, and this sign is important. Without loss of 
generality, we assume t = 1. 

Another symmetry operation of interest is the particle- 
hole transformation Cr i — > (— l)''cj. Under this oper- 
ation, the sign of the nearest-neighbor hopping is un- 
changed while the sign of t' is changed, showing that 
next-nearest-neighbor hopping breaks the particle-hole 
symmetry. Furthermore, the number operator ' — > 
1 — nr as particles are mapped into holes. This symme- 
try was seen in Figures |^ and ^. 

Let us define E(n, t') as the energy per site in the 
ground state of Eq. (^_3|) with next-nearest-neighbor hop- 
ping t' . It is easy to show that as long as n < 1/2, 
E{n,\t'\) < E{n,—\t'\) showing that for small fiUings 
the energy is lower for the positive i'-case. If we in- 
stead consider the region n > 1/2 the particle- hole 
transformation discussed above immediately tells us that 
E{n, —\t'\) < E{n, \t'\), showing that the negative t'-case 
is favorable. At precisely half-filling, the energy is inde- 
pendent of the sign of t' . 

Assume that we consider the case n < 1/2 and t' < 0. 
Then, according to the discussion above, we would gain 
energy if we could somehow change the sign of One 
way of accomplishing this would be to add a uniform flux 
through the system, corresponding to one flux quanta per 
square plaquette. This flux would not affect the nearest- 
neighbor hopping, but it would change the sign of t' and 
therefore lower the energy of the system. Barford and 
Kimo generalized the result of Hasegawa et al. to in- 
clude the t — t' model above and found that in the ther- 
modynamic limit, the kinetic energy is minimized by a 
flux corresponding to one flux quantum per site plus or 
minus one flux quantum per particle. 



VI. HARTREE-FOCK THEORY OF UNIFORM 
PHASES OF Heff 



It has been recognized for some timeE£fE3 that a plau- 
sible response of a Heisenberg antiferromagnet to dop- 
ing, is to form a spiral spin-wave where the Neel order- 
parameter rotates uniformly around a fixed spin-axis as 
one moves along a symmetry axis in the lattice. This to- 
gether with the fact that the doped electron gas favors a 
staggered fiux close to half-filling, indicates that a state 
containing an antiferromagnetic spin-texture that gener- 
ates both a staggered fictitious flux and a spiral-like order 
could lead to an energetically favorable state. 
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In the remaining sections of this paper, we will address 
two related questions. First, we investigate the effective 
Hamiltonian, Eq. (^) , looking for the spin- textures which 
provide the energetically most favorable uniform states. 
In particular, we are interested in whether or not the 
system chooses to incorporate fictitious fiuxes. Further- 
more, it is known that the t — J model has a tendency to 
phase-separate. Concerning the thermodynamic stability 
of the spiral states, Hu et alsn found that in a Hubbard 
model, for small dopings, the spiral phase is unstable 
against phase-separation (for U /t > 10) or domain wall 
formation (for U/t ^ 10). For larger dopings, there are 
regions in the phase-diagram, located around U/t 10, 
where the spiral phases are thermodynamically stable. 
This indicates that, in the t — J model, for small dop- 
ings the spiral state is not thermodynamically stable. An 
interesting question is; if the t ~ J model prefers a flux- 
phase in some region of parameter space, can this ther- 
modynamically stabilize the system, preventing it from 
phase-separation? 

In spir ed by the recent interest in striped phases, in Sec- 
tion VI] we use our approach to model different domain 
walls between Neel-ordered regions. These domain walls 
have an appealing structure as they provide a smooth im- 
plementation of the antiphase boundary and at the same 
time provides the electrons in the doped channel with a 
fictitious flux. 



A. MFT-formulation 

We now look for different uniform phases of the effec- 
tive Hamiltonian given in Eq. (^). The coefficients t^'^ 
and K^^ are now given by 
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-^(1 + cos 6'rr')) 
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where fJr denotes the local staggered spin-orientation 
at site r, and cos^rr' = S^r • f^r'- With this defini- 
tion cos^rr' = 1 for a Neel state and cos 6'rr' — — 1 
for a ferromagnet, where r and r' are nearest neighbors. 
yl(rir, J^r', 2") is the solid angle of the spherical triangle 
spanned by the vectors Or, Or', and z. 

The approach we will use is a simple mean-field theory 
assuming a fixed spin-texture {Or}, defined as the direc- 
tion of the quantization axis, Or. We will assume that 
^rr' — Ox when r and r' are nearest horizontal neighbors, 
and Orr' = Qy when they are nearest vertical neighbors. 
First of all we perform a standard mean-field decompo- 
sition of the Hamiltonian, allowing only for mean-fields 
carrying no charge and momenta zero or Q = (7r,7r). 
This results in the following Hamiltonian: 
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(2 -I- cosOx -t- cos6'j,)A°l - (l-f cos Ox) cos k^/^l^az - [l + cos9y) cosky{Al a2 + Aj? (73] 
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(2 + cos 0. + cos ey){Ay - (1 + cos e,)iAlf - (1 + cos 0,) ((A^,^)^ + (A^^)^) 



(15) 



where we have introduced a two-component vector ^'k = 
(ckjCk+q)', mixing momenta and Q. We have intro- 
duced the Pauli matrices as a basis for the 2x2 matrices 
coupling the ^'k's, although we want to emphasize that 
they have nothing to do with spin in this context. The 
sum over momenta is reduced to half the Brillouin zone, 
defined by BZ' = {\kx\ + \ky\ < n : — tt < kx,ky < 
tt}. Furthewnore, we have only kept those four order- 
parametera23 that turn out to be non-zero numerically. 
These four fields are defined through 

- EkeBZ'(<l*k) 
, ^ly = EkGBZ'COSfcy(*j,a2*k) ^gx 

AL = EkeBZ'Cosfc.(vI/[a3*k) ' 

. ^ly = EkeBZ'COSfcj,(^'j,a3*k) 

where the average ( • ) denotes a thermal expectation 
value with respect to the Fermi-distribution of quasipar- 
ticles of Hmf- The order-parameter A° is simply the 
number of particles in the system, while the other three 
correspond to hopping induced through the term rtr^r' 



in the effective Hamiltonian, Eq. (|6|). In particular, we 
note that A° and A^ are diagonal and hence do not 

cx J y 

mix momenta k and k -I- Q. On the contrary, A^^ does 
mix the two, and therefore carries a momentum Q. As 
can be seen from Eq. (p^), this term only exists in the ki- 
netic term when there is a non-zero staggered flux which 
reduces the translational symmetry of the model. 



B. Instability towards spiral-order at low dopings 

Before turning to the numerical results, let us now dis- 
cuss the electron gas in Eq. (^) at low dopings, confin- 
ing our discussion to coplanar spin configurations and 
neglecting exchange effects in the Heisenberg term so 
that all order-parameters, except A°, in Eq. (|l6|) are 
zero. We choose a spin-structure Or • 0^^^ = cos 6^ and 
Or • O , ^ — COS 9y, allowing for an asymmetry between 
the X- and y-directions. This yields a trivial system which 
is exactly diagonalized. The total energy per site as a 
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function of density n and 0i {i = x,y) is given by 



«k<eF(n) 

-Jn^ (cos6'j; + cosSy + 2) (17) 



1 



where spin) is the Fermi energy corresponding to density 
n and 



sm — cos 
2 



sin — cos ky 
2 



(18) 



Our description of the spin-order in terms of the 9- 
angles does not distinguish between spiral states and 
so called canted states. They both lack fictitious flux 
and both have the same relative angle between nearest- 
neighbor spins. Only second neighbor terms resolve this 
degeneracy. The difference between these two classes of 
states is illustrated in Fig. ^. 
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FIG. 5. The difference between a (11) spiral state (left) 
and a canting state (right) is illustrated. 

If we make a series expansion of the energy per site in 
terms of the doping x we find the following expression to 
first order in x 



Ee^,e,{x) = ~2t (^sin ^ + sin a; - 1 J(l - 2x) 
X (cos 9x + cos 9y + 2) + 0{x'^). 



(19) 



Remember that we have made the transformation 9i i — > 
■n — 9i, expressing the order relative to the antiferromag- 
net instead of the ferromagnet. The energy is minimized 
by 9x = 9y = 2 arcsin[2te/ J(l — 2x)] ~ Atx/J for mod- 
erate dopings. The dependence on tx/J clearly shows 
the competition between the kinetic energy, which drives 
the system towards ferromagnetism, and the Heisenberg 
term, which favors antiferromajgnetism. These results are 
consistent with those of Schulzlil_who also found this in- 
stability. Shraiman and SiggiaEj'Ej, using a more elabo- 
rate method, also found this instability, but their spiral 
state has its pitch-vector along the (lO)-direction rather 
than the (ll)-direction as is found here. In the Hubbard 
model, it is known from Hartree-Fock theory that the 
antiferromagnetic state is stabilized by, the opening of a 
gap at the Fermi surface, and Schula^ has argued that 
a modulation of the spin-order along the (lO)-direction 
opens up a gap more effectively than a spiral along the 
(ll)-direction, i.e. it opens up a gap at a larger part 



of the Fermi surface. The first order theory described 
in this section does not take this fact into account. We 
note that the deviation from Neel-order is proportional 
to the doping and in the limit J 0, the deviation 
from Necl-order becomes large for any pjaiie doping. This 
is consistent with Nagaoka's theorem^, stating that in 
the limit J ^ 0, a single hole doped into a Heisen- 
berg antiferromagnet drives the system into a ferromag- 
netic state. The energy of the spiral state is given by 



E = — j({^2x) ~ "^^^ ~ ^'^'^ hence we have for the 
second derivative of E with respect to the filling n. 



J{1 - 2x) 



3 ■ 



(20) 



We note that the energy versus filling is concave for dop- 
ings X < 1/2, showing that the spiral state has an in- 
stability towards phase-separation for all J > 0. Recall 
that our expansion is only valid for small dopings, and in 
this regime we expect the above conclusion to hold. Also, 
as mentioned earlier, we cannot distinguish between the 
spiral- and canting^states within this approach. Accord- 
ing to Kane et aZo quantum fiuctuations seem to stabi- 
lize the spiral state compared to the canting state. 

A similar analysis can be performed for fillings close 
to zero, i.e. n ^ 1. In this case, measuring 9i with re- 
spect to the ferromagnetic configuration, the analogue of 
Eq. (|9|) becomes (now keeping terms up to second order 
in n) 



Ea 



'2t 



-2Trtn^ \ I cos — cos — 
2 2 

+ ^Jn^ {cos 9^+ cos 9y- 2) + 0{n^). (21) 

In particular we note how the second term introduces 
a coupling between the spin-order in the x- and y- 
directions, allowing for an asymmetry between 9x and 
9y . Minimizing the above energy with respect to the an- 
gles 9i, we find that the ferromagnetic state is stable up 
to a finite doping. At this point different things can hap- 
pen depending on J/t, the system can pick a state where 
9x = and 9y = n (or vice versa), i.e. the system orga- 
nizes itself ferromagnetically along the x-direction while 
being an antiferromagnet along the y-direction. Another 
possibility is that the system chooses a pitch- vector along 
(11) with 9x — 9y — 2 arccos ■ 

The picture we have obtained is therefore that starting 
at zero filling, the system remains in a ferromagnetically 
ordered state up to some threshold value of the filling. 
This threshold increases with decreasing values of the ra- 
tio J /t. Above this threshold filling, the system can be 
a spiral spin-wave with pitch along (10) or (11). When 
the filling approaches n = 1, a (11) spiral state is op- 
timal which continuously merges with the Neel-state as 
n — > 1. But we find in all cases that the system is unsta- 
ble against phase-separation for small dopings. 
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C. Numerics 

Given the thermodynamic instabihty of the spiral spin 
waves, we use the full Hartree-Fock theory and take into 
account uniform phases that have a splay in the Neel 
order-parameter. The aim is to search the space of spin- 
textures, parametrized by {9x,0y,^), to determine the 
one which minimizes the free energy of the system. We 
will consider the following two types of spin-textures: 

1. Coplanar states, described by $ = 0. 

2. TT-flux state, described by 6^ — 9y, ^ = tt. (22) 

The argument for only considering $ = and $ = tt 
states was given in Section ^ Fig. ^ showed that the en- 
ergy per site for a system of free fermions in a staggered 
flux is minimized by either of these two choices. 

Our numerical algorithms work within the grand 
canonical ensemble, with a free energy G{T, fi) = 
(Hmf) — TS — , assuming a fixed chemical potential /z. 
Diagonalizing the mean- field Hamiltonian, Eq. (|l5|), we 
obtain a set of quasiparticle states specified by their mo- 
menta k S BZ', where BZ' is the reduced Brillouin zone 
corresponding to Ifc^; + fcyj < tt, and band-index. The 
band index refers to the two bands occurring because 
of the staggered flux, a S {1,2}. If we then minimize 
the free energy G'(T, fi) with respect to the occupation 
numbers /ka of the quasiparticles, we find that they are 
distributed according to the Fermi-distribution. The en- 
tropy, S, introduced above is defined through, 

S = -ks ^ (^/ka hi /ka + (1 - fka) ln(l - fka)j , 
Q,k6BZ' 

(23) 

where a € {1, 2} labels the two bands. 

In the analysis of the numerical data we would rather 
consider the free energy as a function of the number of 
particles N than the chemical potential /i. This can be 
achieved by forming the Helmholtz free energy through 
the following Legendre transformation; F{N,T) — G + 
^N. Having the Helmholtz energy, we can use the 
Maxwell construction to discuss the thermodynamic sta- 
bility of the phases. 

Analyzing the Hartree-Fock theory involves the follow- 
ing procedure: Given a set of coupling constants, a spin- 
configuration, temperature T, and chemical potential /z; 
pick a set of initial values of the mean-fields. Then solve 
for the quasiparticles of Eq. (^5|) and calculate the new 
mean- fields using Eq. (^). The procedure is then it- 
erated until the mean-fields have converged to a fixed 
point corresponding to a minimum in the free energy. 
Since we are interested in the spin-configuration mini- 
mizing the free energy F{T,N), we will manually vary 
the spin-configuration parameters searching for a global 
minimum of the free energy. 

In practice, rather than choosing a certain value of the 
chemical potential, we choose a fixed filling, successively 



adjusting the chemical potential during the iterations. 
Since the chemical potential was assumed to be fixed 
during the derivation of the self-consistency equations, 
we have to make sure that the algorithm converges to the 
correct fixed point. We have checked explicitly in several 
cases that the correct fixed point is found. An advan- 
tage with this method is that we can access all possible 
fillings, whether or not it is a thermodynamically stable 
region. This is not the case if we specify the chemical 
potential, since the function /x(n) is not invertible. 



7=0.75; r=0.10;L=16 



Co -100 





— Antiferromagnet 




■^-^^^"^ ^ Staggered flux 




Coplanar 







0.2 0.4 0.6 0.8 1 

n 



FIG. 6. The Helmholtz free energy for a,t — J model with 
J — 0.75 at temperature T = 0.1 is shown for the coplanar, 
staggered-flux, and antiferromagnetic spin-configurations. 
Solid lines correspond to Hartree-Fock calculations, while the 
dashed lines correspond to Hartree calculations. 

Let us start by considering a simple Hartree- 
approximation, see Figures |^ and 0. This corresponds 
to putting A^y — A^j. — Aly — 0. Starting at zero tem- 
perature, we find that within the Hartree approximation 
there is a critical value Jc ~ 1 of the coupling J below 
which the coplanar phase always dominates over the fiux- 
generating configurations. When J > Jc, there appears 
a region around n ~ 0.5 where the flux-states are the 
energetically lowest states, see Fig. |^. When the flux- 
state minimizes the free energy, the ^-angles are 7r/2, i.e. 
the maximally allowed values. This state corresponds to 
having the spins distributed along the equator with an 
angle 7r/2 between two successive spins. The solid an- 
gle spanned by this configuration covers half of the unit 
sphere, ensuring $ = tt. Thermodynamically, however, 
it seems to be favorable for the system to phase-separate 
into regions consisting of a hole-free antiferromagnet and 
a hole-rich coplanar structure, respectively. If we con- 
sider smaU dopings {x <g; 1), the 0- angle of the optimal 
state (coplanar) is successively reduced to zero as x ^ 0. 
Since the maximum amplitude of the hopping for a stag- 
gered flux-state is |r| = as determined hy 9 < ^, 
it clearly has a disadvantage compared to the coplanar 
states, supporting |r| = t. The effect of this is that the 
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coplanar state will always be favorable at fillings where 
the kinetic energy is dominant. However, as the filling is 
further increased, the Hcisenberg energy becomes more 
important and it becomes favorable to reduce |r| in or- 
der to improve the Heisenberg energy. At this point, if 
this occurs at a suitable filling, the flux-state can yield 
equally good Heisenberg and kinetic energies, while at 
the same time providing the fermions with a flux that 
can lower the energy eveii further. We know from the 
work of Hasegawa et alu that the energy of a system 
of free electrons on a square lattice is minimized when 
there is exactly one flux quantum per particle. Since our 
flux-state carries a flux tt per plaquette, it will be most 
suitable close to n — 1/2. The coplanar state on the 
other hand carries no flux, and will therefore be optimal 
when n = or n = 1. This competition explains why 
the coplanar state becomes energetically favorable again 
when moving from a flux-phase towards higher flUings. 



7=2.00; r=0.10;L=16 
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FIG. 7. The Helmholtz free energy for a,t — J model with 
J — 2.00 at temperature T — 0.1 is shown for the coplanar, 
staggered-flux, and antiferromagnetic spin-configurations. 
Solid lines correspond to Hartree-Fock calculations, while the 
dashed lines correspond to Hartree calculations. 

In our calculations, a slight technical point should 
be mentioned. When the system phase-separates into 
a coplanar part and an undoped antiferromagnet, the 
chemical potentials in the two subsystems are not equal. 
This can be seen, for instance in Fig. ^, where the 
Maxwell construction connects the antiferromagnetic 
TT-i = 1 point with a point uq ~ 0.55 in the coplanar 
phase. Since ni is an end-point of the free energy curve, 
the derivatives at uq and ni arc unequal, and hence the 
chemical potential is different in the two phases. 



D. Results from Hartree-Fock calculations 

We now apply the Hartree-Fock scheme to the problem 
to understand whether or not exchange effects can resolve 



the near degeneracies found in the Hartree calculations. 
The spin-textures considered are those given by Eq. ([2^). 
Before starting the full calculations, we determined nu- 
merically which of the different order-parameters are the 
important ones. Initially there were tenEj but numer- 
ically we find that only four are non-zero in the spin- 
configurations we have examined. To speed up the calcu- 
lations we have set the other six to zero by hand. We find 
that besides the filling (A°), A^^, and A^^ are important 
for the fiux-phases, while A^^ and A^^ are important for 
the coplanar configurations. 

Numerically we find that the 7r-flux-phase converges to 
a state where the non-zero order parameters are A^y = 
= A, which leads to a quasiparticle spectrum of 
the form Ck = — /i ± 2t-\/cos2 kx + cos^ ky, where fi and 
t are renormalized values of the chemical potential and 
hopping amplitude, respectively. The momentum k be- 
longs to the reduced Brillouin zone. At half-flUing, this 
dispersion relation has four gapless Dirac-points where 
the energy vanishes linearly. These points are located 
at {kx,ky) — (±7r/2, ±7r/2). Similarly, for the coplanar 
states, the mean-flelds renormalize the hopping ampli- 
tudes so that Ck = — A ± cos cos Ky). 

If we look at a typical free energy plot, such as Fig. ||, 
we find, as in the Hartree case, that for low fillings the 
coplanar configuration is the optimal, where at low dop- 
ings it merges with the antiferromagnet. For small J's 
the coplanar configuration clearly dominates over the 
flux-conflguration for all fillings up to the point where 
they merge with the pure antiferromagnet. Increasing 
J brings the fiux-configuration energetically closer to 
the coplanar configuration. However, in contrast to the 
Hartree-case, it does not seem like the flux-state will 
become energetically favorable over the coplanar-states. 
Concerning phase-separation, the picture is very much 
the same as the one described above. For low tempera- 
tures and dopings smaller than roughly 0.5, the system 
favors a phase-separation into parts consisting of a hole- 
free antiferromagnet and a hole-rich coplanar state with 
doping X ~ 0.5. 

Figures ^ and |^ show the free energies for the best 
coplanar- and flux-conflgurations plotted together with 
the free energy of the pure antiferromagnet in two differ- 
ent cases, J — 0.75 and J = 2.00. The energy-scale is 
flxed hy t — 1. The temperature is set to T = 0.1 and 
the size of the system being considered is 16 x 16 sites. 

In Fig. 0, it is clearly seen how the flux-phase dom- 
inates over the coplanar-phase close to n ~ 0.5 at the 
Hartree-level, but not in the Hartree-Fock approxima- 
tion. 

We have generalized the t — J-model by including 
a nearest-neighbor Coulomb repulsion through a term 
r' '^r'Tr'- In our effective model this corresponds to 
redefining i^rr' ' — ^ K^r' + V . When including this term 
the order-parameter corresponding to a charge-density- 
wave, A^ = X]keBZ'(^k'^i^k), becomes important. As 
it turns out, a positive value of V can favor the flux- 
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phase as is seen in Fig. |^ where we have shown the free 
energy versus filhng for a system described by J = 1.25, 
K = W/J = 2.00, and temperature T = 0.1. As can 
be seen from this figure, there appears a narrow region 
around n ~ 0.53 where the flux phase has the best ener- 
getics. 

7=1.25; A:=2.00; r=0.10; L=16 




FIG. 8. The Helmholtz free energy for a t — J 
model with J = 1.25 and K — 2.00 at temperature 
T = 0.1 is shown for the coplanar and staggered- flux 
spin-configurations. Note the small region about n ~ 0.53 
where the fiux-configuration is energetically more favorable 
than the coplanar spin-configuration. 
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FIG. 9. This figure shows the state which minimizes the 
free energy depending on doping and the coupling constant 
K. We have chosen J — 1.25, temperature T = 0.10, and 
the finear dimension of the system is L — 12. The data is 
based on Hartree-Fock calculations. As can be seen, a narrow 
region where the flux-configuration is the best uniform state 
appears around n ~ 0.53. The coplanar phase continuously 
merges with the antiferromagnet as n ^ 1. 



Finally, in Fig. ^ we show a phase-diagram as a func- 
tion of filling and coupling K. In this data we have 
fixed the exchange coupling J = 1.25 and temperature 
T = 0.1. The figure shows how a narrow region of a 
staggered- flux phase occurs close to n ~ 0.5. 

Summarizing our numerical results, we find that a 
staggered flux-phase can be energetically favorable com- 
pared to a coplanar spiral state. However, in the region of 
doping where this happens, the system generically seems 
unstable against phase-separation. Even if this indicates 
that the flux-phase does not occur as a thermodynam- 
ically stable phase in our effective model, it has good 
energetics and it would be interesting to investigate the 
effects of spin-fluctuations in this picture. 

In the next section we discuss stripes within this frame- 
work. We know that stripes are antiphase boundaries be- 
tween antiferromagnetically ordered regions. One way to 
model this is to have a spin-order twisting as the bound- 
ary is crossed. Since we also know that the stripes are 
doped, it is tempting to think that the twisting is such 
that a flctitious flux is generated. The picture is also ap- 
pealing since it makes use of the instability of the doped 
antiferromagnet towards phase-separation. We will re- 
turn to this topic in Section VII 



E. Circulating currents 

An important issue to address is whether or not the 
flux-states are accompanied by circulating currents in the 
system. To answer this question, we consider the current 
operator at site r in the (5-direction, ^^(r), which can 
be identifled from charge conservation together with the 
Heisenberg equation of motion. 



[jsir)-jsir~6)] 



dp{r) 
dt 



The result is a current operator taking the form 



h 



h.c. 



(24) 



(25) 



We decompose the current into uniform- and staggered 
parts as js{j^) = Jsi'"^) + Using the mean-field 

decomposition we find that the expectation value of the 
uniform currents vanish. The uniform currents are pro- 
portional to some of the order-parameters that have been 
left out of the discussion. We know that our Hamiltonian 
is invariant under a lattice translation followed by a time- 
reversal operation, but this composite operation reverses 
the direction of the uniform currents which hence must 
vanish in a thermodynamic expectation value. The stag- 
gered currents take the form 



(J^(r)) 



2t 



2 ' 



10 



t),'(r» 



(28) 



The currents are gauge invariant, and the formal lack of 
symmetry between the currents in the x- and y-directions 
is due to gauge choice. The order-parameters also are 
gauge dependent, and this restores the symmetry be- 
tween currents in the two directions, see Fig. QG. 
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FIG. 10. A site in the even sublattice is shown. Current 
conservation takes the form 2{jx + jy) = 0, which we also 
observe in our numerics. 



According to the previous subsection the non-zero 
order-parameters in the 7r-flux state are Aj?^ and A^^, 
and we conclude that there are no circulating currents 
in the tt-Aux state. Obviously, the flux-free case is also 
free of currents, as can be seen from the equations above, 
in this case it is A^^, rather than A^y which is non-zero. 
Noting that in the case $ = ±7r the Hamiltonian is time- 
reversal invariant so the the one-particle states possibly 
carrying current must be degenerate in energy. In a ther- 
modynamic ensemble, all states of the same energy are 
weighted equally, and the physical currents in a system 
cancel and there will be no staggered currents in the sys- 
tem. For any flux not being an integer multiple of tt, this 
symmetry is lost and circulating currents appear. 

If we calculate the above expectation values in a state 
having a flux <I> ^ mr, we find circulating currents in the 
system. 

A proper treatment of this problem should consist of 
a gauge invariant coupling to a real electromagnetic flux 
in addition to the spin-generated fluxE3. Let ^mag de- 
note the real electromagnetic flux through a plaquette, 
induced by circulating currents. The energy cost of cre- 
ating the magnetic fleld should be added through a term 



J^mag — „ 



2 ' mag ■ 

plaquettes 



(27) 



The constant K is given hy K — ^J^^i^i , where a is the 
two dimensional lattice constant, d the distance between 
the copper oxide planes, and [i ^ th e relative permeabil- 
ity. As was discussed in Section III , this constant is huge 
compared to the typical electronic energies. The total 
energy as a function of the electromagnetic flux is then 
written as 



where -Et- j($mag) denotes the energy of the t — J model 
when there is an extra flux <&mag in addition to the spin- 
generated flux. Minimizing the energy with respect to 
the electromagnetic flux leads to an equation of the form 

.g(^'mag) + Ji^^mag = 0, whcrC g(*mag) = £^t'_ j(*mag)- 

As we have pointed out, the magnetic energy-scale K is 
much larger than the electronic energy which is of the or- 
der max(t, J). As a consequence of this, the magnetic flux 
will be suppressed and it is reasonable to put <i>inag = 
when we solve for the eigenstates of the system. For this 
to be a self-consistent solution there must be no currents 
in the system due to the spin-generated flux. This is true 
for the cases $ = or tt which we have focused on. 



VII. STRIPED STRUCTURES 

Striped structures forming antiphase domain walls be- 
tween undoped antiferromagnetic regions have been ex- 
perimentally observed in the doped high-temperature 
superconductorsE3. We would like to understand if a 
striped phase can be explored using the effective t — J 
Hamiltonian, Eq. (^). There are several facts that make 
this an appealing approach. We have already seen in 
the previous section that for low dopings, the uniform 
(spiral) states are unstable against phase-separation. Us- 
ing our spin-polarized approach we can create a smooth 
antiphase boundary, successively changing the order- 
parameter from to —z. And, since all the holes of the 
stripe are located in the domain wall it could be favor- 
able for the system to generate a flux in this region. This 
can be accomplished using a spin-texture, which simul- 
taneously generates the antiphase boundary. Further- 
more, the experimentally observed value of the doping of 
stripes injI«|a2_a;Sra;Cu04 is 0.5 holes per unit length of 
the stripet^. This is close to the region where we have 
seen that a 7r-flux may be favorable. (See Section 0.) 
Inspired by these nice properties we have undertaken an 
investigation of striped phases within our approach. Our 
main ambition has been to gain an understanding of what 
such a striped phase would look like, and, in particular, 
whether a fictitious flux is exploited or not in our model. 
The technique we use is a self-consistent Hartree calcu- 
lation. 
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FIG. 11. Two possible spin-orders making up an antiphase 
boundary between two Neel-ordered regions. The y coordi- 
nate is orthogonal to the direction of the stripe, and y = G 
corresponds to the center of the stripe. The left part of the 
figure shows the antiferromagnetic order-parameter passing 
through zero, while the right part shows our scenario with 
the order-parameter tilting as one passes through the stripe. 
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Let us consider a system with an antiphase boundary 
along the x-axis, and with the antiferromagnetic order 
being +'z at y = oo and ~z at y = —oo. Two possible 
scenarios for an antiphase spin-order in a stripe comes 
to mind. First, the amplitude of the antiferromagnetic 
order-parameter may simply decrease, passing through 
zero and becoming negative as one passes through the 
stripe. This scenario preserves the rotational symmetry 
about the spin 2-axis. The second possibility is that the 
Neel order-parameter starts to tilt as one approaches the 
stripe along the y-axis, lying within the spin x-y-plane 
at the center of the stripe and then rotating towards the 
positive z-axis as one goes to y = cx). Introducing holes 
in the domain wall, the amplitude of the spin-order will 
decrease as it depends on the particle density through 
Sr = ■^rir^lr, but it wiU not vanish. An illustration of 
these two scenarios is shown in Fig. |ll|. From an ex- 
perimental point of view, Tranquada and coworkers ar- 
guerlri that there are indications speaking in favor of the 
first scenario, but these indications are not conclusive. 
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FIG. 12. To the left is shown a site-centered stripe, 
in which there is a site at which the spin-order is in the 
X- J/- plane. In the bond-centered case (right), the antiferro- 
magnetic spin-order resides in the x-y-plane at an imagined 
point between two sites. 

Four different stripe geometries with preserved uniax- 
ial symmetry come to mind. First, stripes can go along 
either the (10)- or (ll)-directions and we can choose ei- 
ther site- or bond-centered stripes, all in all four possible 
combinations. The difference between site- and bond- 
centered stripes is shown in Fig. 
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FIG. 13. Our model of the striped phase is shown. A and 
B denote the two antiferromagnetic phases (order-parameter 
±1) and in between each A and B pair is an antiphase 
domain wall reversing the sign of the antiferromagnetic or- 
der-parameter. The figure also introduces ^ as the width of 
the antiphase boundary, and d as the size of the antiferromag- 
netic layers between the stripes. 

Our aim is to find the distribution of holes and spin- 
texture that minimizes the energy of an antiferromagnet, 
when we assume the holes are arranged in stripes. The 



configurations we consider have the structure shown in 
Fig. |l3|, consisting of a repeated structure of Neel-regions 
and antiphase domain walls. The order-parameter of the 
antiferromagnetic regions change sign every time a do- 
main wall is passed. 

We will assume that all the holes are located near the 
domain walls. Following the notation in Fig. we will 
denote the width of a single antiphase boundary by ^, 
and the width of an antiferromagnetic region by d. The 
average particle density in the domain wall will be de- 
noted Ud, and we will also use the number of holes per 
unit length (along x) of the domain wall, (5 = (1 — nd)£,- 
Assuming an average filling n of the system, we have 
{i + d)n ~ + d. 

Before writing down the total energy of this configura- 
tion, we make the following observations. We note that 
the system in Fig. ^ is build up from units (antiferro- 
magnets and domain walls) and we would like to write 
the total energy in terms of the energies of the individ- 
ual units. To do this we define the energy per site of 
an antiferromagnetic unit, Eaf, as the energy per site of 
the antiferromagnet with periodic boundary conditions 
in the x- and y-directions. The energy per site of the do- 
main wall as a function of the number of holes per unit 
length of the stripe, Ed{S), is similarly defined by putting 
periodic boundary conditions on the domain wall. In this 
case it is important that the edge of the domain wall is 
Neel-ordered and undoped. If this condition is not ful- 
filled, there will be surface energies associated with the 
gluing of a domain wall to an antiferromagnetic region. It 
is easy to see that minimizing the total energy of the sys- 
tem is equivalent to minimize the energy per introduced 
hole in the domain wall. 

Our approach will therefore be to consider a single do- 
main wall and minimize the energy per introduced hole 
with respect to the parameters describing the domain 
wall. These parameters contain spin-texture related pa- 
rameters (which will be introduced shortly) , the number 
of holes per unit length of the stripe (5), and finally we 
have the four options for the stripe geometry; site/bond- 
centered and direction (10)/(11). 

Note that the definition of ^ is somewhat arbitrary 
in the sense that there in practice may be a smooth 
crossover from the domain wall to the Neel-ordered re- 
gion. For this reason, (5 is a better measure of the stripe- 
doping than Ud which depends on ^. Knowing S and ^, 
the stripe periodicity is 



l = 2iC + d) 



2S 

X 



(29) 



where we have introduced the average doping of the sys- 
tem, X — 1 — n. Thus we find that for low dopings, the 
separation between the stripes scales as x'^ . This de- 
scription is valid as long as the stripes remain separated 
so that we can neglect stripe-stripe interactions. This 
condition is fulfilled as long as x <C |. We note that as 
we change the overall doping of the system, the structure 
of the isolated stripes remains, at least as long as d ^ 1. 
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Thus the wavevector describing the spin order is 27r/Z, 
and the wavevector of the charge order is twice that, i.e. 
At: /I. 



A. Stripes in the (lO)-direction 

Let us start with a description of how a single (10)- 
stripe is modeled. Although we are interested in infinite 
domains, in the numerical simulations we are forced to 
work with finite stripe width, which we denote by the 
integer w. We will assume the local spin orientation of 
the stripe to described by a unit-vector field fir as 

51r = (— l)'"(sin 6'r cos 0r, sin 0r sin 0r, cos ^r), (30) 

where the spherical angles Or and 4>r are functions of po- 
sition. The spatial dependence of these angles is assumed 
to take the form 



arccos 
qxX. 



tanh ■ 



l-2y 



2C 



(31) 



The construction can be thought of as a cylindrical pro- 
jection (also known as Mercator projection) of the lattice 
into the spin sphere by associating latitude and longitude 
with the X- and ^/-coordinates, respectively. 

L^=2%lq^ 



t t t t t t t t t t 
t t t t t t t t t t t t 



,.i....J....i....i....i....i...J....i....i....i....i....i. 



\ \ \ \ \ 



\ \ \ 



1 1 [ i i i [ i i [ [ L 



FIG. 14. A sketch of the Neel order-parameter in the 
neighborhood of a stripe constituting an antiphase boundary 
between two regions with perfect Neel order. All spins are 
unit vectors. In the figure we have also indicated the length 
scale 5. Lx, the period along the x-direction, is related to 
qx through = 2'k/Lx- This particular spin-configuration is 
clearly site-centered. 



We are interested in the limit w — » 00. This limit has 
two discrete values, depending on whether w is even or 
odd. For odd values of w the spin-order of Eq. ( ^l| ) is such 
that the spins at y = (ii;-|-l)/2 are lying in the a;-?/-plane, 
and hence the stripe is site-centered. If we instead con- 
sider an even integer w, the center of the stripe is located 
in between two rows and the stripe is bond-centered. 

In addition to the discrete choice of bond- or site- 
centered stripe, there are two free continuous parameters 
in Eq. (|3l|); ^, which determines the characteristic width 
of the spin-domain wall, and qx, which is the pitch of the 
rotation of the spin about the z-axis along the length of 
the stripe. Fig. |lj illustrates this construction. 

Making an apparently trivial point, which we return 
to in the next section, we note that there is no choice 
for ^ and qx which unwraps the antiferromagnet into a 
ferromagnet; taking = tt will create a ferromagnetic 
channel in the center of the domain wall, running along 
the stripe, but will not affect the z component of the Neel 
order-parameter . 

When writing down the Hamiltonian of this system we 
need to fix a gauge. If the spin-configuration around a 
square plaquette in terms of spherical angles is given by: 
{9i,(t)i), (7r-6'2,<^i-H7r), (6*2, (/'2), and (7r-6li, 02 +7r), the 
product of the hopping elements around the plaquette is 



T'0lTi2T23'''30 



- sm t'l sm t'2 sm 



2 "^1 



72 .2 91 

— sm 



<0, (32) 



showing that the flux through a plaquette is exactly equal 
to TT. The texture is therefore quite natural for the fol- 
lowing reasons: It is periodic and has a uniform flux per 
square. The flux tt per plaquette favors heavily doped re- 
gions near the center of the stripe, while the effective hop- 
ping in the Neel regions, where the system is undoped, 
is vanishing. With flux ±7r through the plaquettes there 
is no broken time-reversal symmetry, and there are no 
circulating currents or induced local magnetic fields. 

In order to perform our calculations, we use a mixed 
representation using momentum space in the x-direction 
and real space in the y-direction due to the transla- 
tional invariance of the system along the cc-direction. The 
length of our system along the stripe will be denoted by 
L. Furthermore we assume that the number density at a 
site r only depends on the y-coordinate, i.e. n(r) = n{y). 
We will use the following Hartree-decomposition of the 
interaction term in Eq. (|6|) 



TirTir' ~ n{y)nr' + iivn{y') - n{y)n{y'). 



The effective Hamiltonian can then be written as 



w—l 



(33) 



i?buik = -2t'Y^[~iy^cos{k)Tx{y)nk,y ~t^Ty{y)^ [4,^+1^;.,^ + h.^ 

y—l k y—1 k 

J 

+ jY.(hx{y) - 1) 51 ['^nk,vn{y) - n{yf] 



v=i 
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J 



7 (^y(y^ ^ 1) + 1) + '^fc.a+i'^Cy) - n{y)n{y + 1)], 



(34) 



where we have introduced the coordinate-dependent cou- 
pUng constants 



(3,5) 



Tx{y) = 



(35) 



and the analogous relations for hy{y) and Ty{y). The 
factor (—1)^ in the hopping term of Eq. (Q) comes from 
the flux ±7r through each plaquette. In addition to the 
terms in Eq. ( ^^ , we have also added a local potential at 
the upper and lower edges of the system to simulate the 
effect of the adjoining antiferromagnetic region. Without 
this potential, which has the form 



J 



(36) 



the system may gain energy by expelling the holes to 
the edges where they break fewer antiferromagnetic links. 
The total Hamiltonian is then given by 



bulk 



(37) 



The numerical calculation involves solving self- 
consistently for a charge profile described by n(y), where 



(38) 



and the expectation value is with respect to the Fermi- 
distribution of quasiparticles of Hiq. We use a chemical 
potential to control the overall number of particles in 
the system. The chemical potential is determined during 
each iteration of the self-consistency equation, Eq. (38). 
Our calculations are performed at temperatures close to 
zero, T = O.Olt, and the quantity we focus on is the free 
energy as a function of the number of particles and the 
parameters describing the stripe, F{T, N) = (i?io} — TS. 



B. Stripes in the (ll)-direction 

The analysis of the (ll)-stripes is similar, although the 
geometry of the stripe introduces some complications. In 
Fig. |l^, we have shown the geometry used for these con- 
figurations. Due to the tilting of the lattice, we note that 
every second row of constant y in the lattice is shifted 
to the right by half the lattice spacing in the x-direction. 
Moreover, if we denote the even- and odd sublattices by 
A and B, respectively, we find that all points belonging 
to A reside on points having odd values of y, while those 
belonging to B are assigned even values of y. This is 
shown explicitly in Fig. hq. 




.(4,4) 



,(4,2) 



The geometry used for the description of 
(ll)-stripes is shown. The original square lattice has been 
tilted 45 degrees, and we have labeled each point with a co- 
ordinate {x,y). These coordinates are not be confused with 
the coordinates in the original, untilted square lattice. Each 
lattice point has also been marked with the sublattice, A or 
B, to which it belongs. The domain wall is directed along the 
X-axis in this coordinate system, and coupling constants will 
depend on the y-coordinate only. 

The stripe is directed along the x-axis in the coordi- 
nate system defined in Fig. |l^, and the spin-configuration 
is given by 

rir — (— I)^(sin6'rcos0r,sin0rsin(/)r,cos0r)- (39) 

We note that 0^ = Q and 0^ = 1^ correspond to the two 
Neel-states, and hence we can (as in the (lO)-case) in- 
terpolate between the two by continuously changing 9^. 
across the domain wall. To be explicit, we will use the 
following parametrization of the spherical angles: 



tanh ■ 



l-2y 



2C 

q^x y odd 

qx{x + 1/2) y even 



(40) 



Note that we have shifted x by 1/2 for even y- values to 
account for the shift of lattice points in the x-direction 
as discussed above. It is also important to stress that, 
contrary to the (lO)-case, it is possible to recover a fer- 
romagnetic configuration by a suitable choice of param- 
eters, namely making ^ large and taking q,j. = lix. This 
corresponds to Q-r = 7r/2 and q^ — 27r for odd j/'s and 
qx — TT for even y's, respectively. The effect of the rota- 
tion due to q^ is therefore to rotate all spins belonging to 
sublattice _B by tt about the spin z-axis. The result is a 
ferromagnetic configuration, where all spins point along 
the positive x-axis. This difference between the (10)- and 
(ll)-stripes reflects the fact that the local field along ei- 
ther side of a (ll)-stripe is ferromagnetic whereas in the 
(10) case it is antiferromagnetictZl. 
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The next issue we will address is the properties of the 
fictitious fluxes generated by a certain spin-configuration 
in the domain waU. Fig. |l^ shows the flux-pattern which 
is generated from the spin-configuration in Eq. (HQ). 




FIG. 16. The lattice of the (11) domain wall is shown. 
The structure of the flux pattern is given by the $y's. In 
particular we find that the flux is uniform along the x-axis. 
We have fixed a gauge by defining the phases associated with 
the links in the lattice. The gauge is chosen such that all 
non-zero phases are on links connecting sites with the same 
a;-coordinate. 

We find that the flux is uniform along the x-direction, 
while it depends on the y-coordinate of the plaquette to 
which it belongs. Fig. |l6| also fixes a gauge, defined by 
the phases (j>y. As we can see from the figure, there is one 

Hbulk - E E {ry{y)e'*^'^ + Trf(y)e^(-i) 

y=l k 
w — l J 

+ E 4 C^v^y^ + ''diy) - 2) E {r^k,v 
y=i k 



phase more than there are gauge invariant fluxes. There- 
fore, we can choose a gauge in which (pi = 0. Next, we 
consider the flux 4>i which is determined by $1 — (f>i+<j)2- 
More generally, we find = {—l)^^-^ {(py + (f)y-^-l) . In this 
way we can solve for the phases 4>y in terms of the fluxes, 
finding 

cby+i = {-i)y E (41) 

y' = l 

Numerically, given the spin-configuration of the domain 
wall we calculate the fictitious fluxes, $y, using Eq. (|^) 
and then use Eq. (|l]) to find the appropriate phases that 
enter the Hamiltonian. 

We define the amplitudes of the hopping and Hciscn- 
berg interactions as 

hd{y) — ^x,y ■ ^x+(-l)y,y+l 

which allows us to write the bulk part of the Hamiltonian 
as 

T'iy + 1) + nk,y+in{y) - n{y)n{y + 1)] . (43) 



For the same reasons as in the (lO)-case, we will add a 
local potential to the vertical boundaries (note that each 
boundary site connects to two sites in the environment) 

fe 

The total Hamiltonian, Hn, is the sum of the bulk- and 
boundary-contributions, i.e. 

Hii = -ffbulk + -H^b.c- (45) 

The numerical procedures are completely equivalent to 
those used in the (lO)-case. 

C. The optimal stripe 

According to our model, the physically relevant stripe 
configuration is that which minimizes the domain wall 
energy per introduced hole. We will use the undoped 
antiferromagnet as an energy reference state, this being 



the optimal state at zero doping. As we dope holes into 
the system, the total energy of the domain wall will be a 
discrete function of the geometry (site- or bond-centered 
and direction (10) or (11)) and a continuous function of 
Nh, cjx, and ^. Note that since we work with a chemical 
potential, the number of holes in the domain wall, Nh, is 
not restricted to be an integer. 

We label this energy E{Nh, qx,£,)- However, the physi- 
cally interesting quantity is the number of holes per unit 
stripe length, S — Nh/L. We define the domain wall 
energy per hole according to 

Eh[o,qx,t.)^ ■ (46) 



As we argued in the beginning of this section, to find 
the optimal stripe-configuration we have to minimize this 
function with respect to 6, Qx, and ^ for bond- and site- 
centered (10)- and (ll)-stripes. 
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indicate that J ~ 0. We also find that ^ is very small 
for this optimal stripe, i.e. there is a sharp spin-domain 
wall with a single tilted row of spins. This row is ferro- 
magnetically ordered as we find = vr, and the holes 
are tightly confined in the neighborhood of the domain 
wall. The spin- and charge-profiles are shown in Fig. 

Since the optimal domain wall is so narrow, the prod- 
uct of the effective hopping-amplitudes around any pla- 
quette in the lattice is approximately zero. Thus we must 
conclude that the system does not take advantage of the 
fictitious TT-fiux through the plaquettcs. 

As Fig. |l^ shows, the bond-centered (ll)-stripe is en- 
ergetically very close to the optimal (lO)-stripe described 
above. An illustration of this domain wall configuration 
is shown in Fig. This (ll)-stripe is characterized by 
6 ~ 0.71, Qx ~ 2tt, and f ~ 0.73. It is important to 
point out that qx — Stt is not equivalent to Qx = since 
Qx — 27r performs a 7r-rotation about the z-axis of one 
of the sublattices, as we discussed in the previous sub- 
section. Furthermore, we want to emphasize that this 
stripe-configuration does not induce any fictitious fiuxes 
Turning to our numerical results, we have initially con- ^"^^ consequently there are no currents that could ener- 
sidered a system with a Heisenberg coupling J = 0.40, gctically disfavor this configuration, 
where the energy scale is fixed by i = 1. This value 

was chosen because it corresponds to a value of the ex- /-^ X ,X X X X Undoped 
change coupling constant which has been used by oth- "^-^ ''-j^'' q 5% doping 

ers to model the high-temperature superconductors. In ,-' ,-' ,-' , ' \ ,-' ,- 

Fig. ^ we show the optimal energy per hole as a function O 31% doping 

of (5 for the four stripe geometries, i.e. we have plotted (3^ gf jS' ,!Sr £2^ O 31% doping 

En{b)^Y^v^^^,i,Eu{&,qx,i). W W W o 5o/„ doping 
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I I I I I I ' FIG. 19. The structure of the optimal confined (ll)-stripe 
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I T I T I T I H y spins, and the radius of the circles the amount of hole-doping. 
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-^-^--^-^--^-^--^O 36%dopi ng correspond to 31 % hole-doping. Regions being approximately 
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FIG. 17. The energy per hole is shown as a function of 
the number of holes per unit stripe length, 5. An investi- 
gation of the curves shows that the minimum occurs for the 
site-centered (lO)-stripe at 5 ~ 0.46. 
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FIG. 18. The structure of the optimal (lO)-stripe for 
J = 0.40 is shown. Arrows indicate the polarization of the 
spins, and the radius of the circles the amount of hole-doping. 
Small circles correspond to 4 % hole-doping, while large cir- 
cles correspond to 36 % hole-doping. Undoped regions lack 
circles. 

From Fig. ^ we can read off the optimal stripe con- 
figuration, which will make up the domain walls in the 
striped phase. As is indicated in this figure, the opti- 
mal domain wall is a site-centered stripe along the (10)- 
direction having 6 ~ 0.46 holes per unit length of the 
stripe. This agrees with result|S-|from DMRG calcula- 
tions by White and collaboratorg23, who find stripes with 
(5 = 0.5 for J = 0.35. Furthermore, the experimental data 



Since the optimal (lO)-stripe is very close in energy 
to the (ll)-stripes it is interesting to investigate what 
happens as we tune the strength of the Heisenberg inter- 
action, J. Numerically we find that slightly increasing 
J above J — 0.40 favors the (lO)-stripes compared to 
(ll)-stripes, while decreasing J favors the (ll)-stripes. 
There will be a crossing point slightly below J = 0.40, 
at approximately J ~ 0.36, where the (ll)-stripes have 
lower domain wall energy than the (lO)-stripes. 

There is a technical point which is important when 
looking at the energy per hole as a function of (5 for the 
(ll)-configurations. As we mentioned in the previous 
subsection, it is possible to unwrap the (ll)-domain walls 
into ferromagnets. If we follow the (ll)-curves in Fig. |l^ 
to larger values of b we find that the energy successively 
decreases below what we called the optimal stripe con- 
figuration. The configurations that correspond to these 
low-energy states are close to ferromagnetic and with an 
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almost uniform charge-distribution. Physically this cor- 
responds to a global phase-separation into an undoped 
Neel region and a heavily doped ferromagnetic region, 
i.e. the phase-separation discussed in Section VI. Let 



us therefore compare the energetics of the optimal (10)- 
domain wall and the global phase-separation. 

At the minimum of Eh{5) in Fig. |l^, we can read 
off the energy per hole of the optimal domain wall as 
i?dw — —0.76. This is to be compared to the energy per 
hole for the totally phase-separated state E-ps — —1.06 
for the same value of J. This number was obtained 
from a Hartree calculation using Eq. ( |l7|) on a system 
of the same size as the stripe grid. This clearly shows 
that, within our approximation at least, phase-separation 
is energetically advantageous compared to domain wall 
formation. In a real system, the energy of the phase- 
separated state is raised due to the Coulomb interaction 
between the holes and this could make the domain wall 
thermodynamically stable. 

The cuprates seem to favor the formation of rather 
narrow stripes, not phase-separation. We will take the 
point of view that there is some mechanism, not captured 
in our approach, such as a long-range Coulomb interac- 
tion, which prevents grouping all the holes together and 
instead favors the formation of stripes on some interme- 
diate length scale. Therefore, we will only consider the 
stripe-configurations which are local in nature. 



system has t' > 0. In an antiferromagnet, the presence 
of a second-neighbor hopping is important since it allows 
for holes moving through the sublattices without disrupt- 
ing the antiferromagnetic order. Typically, the value of 
t' used in the literature for describing a hole-doped anti- 
ferromagnet is t' ~ —0.3. 

If we consider the limit in which the nearest-neighbor 
hopping is completely frozen out, and there is only 
second-neighbor hopping, i.e. 

Ht.-J = -t' [Tr.'clcr+\i.' 
((rr'» 

+ -^^{^r-?i-r' -l^nj-Ur', (48) 
(rr') 

we note that the transformation i — > (— l)^Cr leaves 
the Heisenberg term unchanged, while the second- 
neighbor term changes sign. This shows that in this limit 
the asymmetry between ±t' vanishes. Hence, the sign of 
t' is irrelevant in the Neel-ordered regions, and it is only 
in the region where the spin-twist occurs that t and t' are 
simultaneously present and accordingly, the sign of t' is 
important. 

Introducing next-nearest-neighbor hoppings as in Sec- 
tion [V allows for new closed particle-orbits in the lattice. 



and hence also for new gauge invariant fluxes. There are 
four of these fluxes and they are defined in the right inset 
of Fig. |0|. 



D. Including next-nearest-neighbor hopping 

An unphysical feature of our simulation is that since 
hopping between the antiferromagnetic spins is forbid- 
den, the effect of second-neighbor (diagonal) hopping be- 
comes important in the Neel state. In our approxima- 
tion of the t — J model, hopping is frozen out for this 
state, so that electron transport will be dominated by 
any second-neighbor terms if they are non-zero. This 
term will permit electrons to diffuse into the Neel state 
and could therefore be expected to delocalize the holes 
from the stripe center. We will extend our model of the 
(lO)-stripes by adding this hopping to the Hamiltonian 
through a term 



-t' 



E 

((rr')) 



'rr'^r' 



Cr -t- h.C. 



(47) 



where ((rr'}) denotes next-nearest-neighbor pairs. As be- 
fore, the calculations reported below are performed with 
the Heisenberg coupling J — 0.40 and at temperature 
T — 0.01 in terms of energy units set by i = 1. In this sec- 
tion, we will only consider the effect of second-neighbor 
hopping on the (lO)-stripe, since this was found to be the 
optimal stripe-configuration for J = 0.40. 

It has been argued in the literature that the sign of 
t' in the high- Tc cuprates depends ctfi|-si/hether the sys- 
tem is hole-doped or electron-dopedE§Ej. A hole-doped 
system corresponds to < 0, while an electron-doped 




FIG. 20. Our gauge choice for the next-nearest-neighbor 
hopping. The gauge choice on the horizontal- and verti- 
cal links is the same as the one used previously for the 
(lO)-domain walls. The right inset defines the fiuxes <b\ 
through the four sub-plaquettes of a square plaquette. 

Investigating the spin-structure of Eqs. ( ^0| ) and (|3l]), 
we find that the four fluxes associated with row 

y are determined by a single parameter '^y according to 

(49) 
wc find 

(50) 



Recalling our topological constraint, Eq. (| 
that it is satisfied as 



27rn, 



with ri = ±1, which is what we expect for an antiferro- 
magnet. Which sign that applies to a certain plaquette 
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depends on the sublattice associated with the plaquette 
as well as the spin-configuration. Fig. 20 defines a gauge 
by introducing the phases ay and fiy. Using Eq. ( |4g| ) it is 
easy to read off the parameters ay and [3y from Fig. |2^. 
Doing this we find ay = —f3y = ^'y for odd values of y, 



and a„ 



-j3y — ±TT — 'i'y for even j/'s. 



As in the nearest-neighbor case, we will work in mo- 
mentum space in the x-direction, while keeping the real 
space description in the j/-direction. We note from Fig. ^ 
that the phases of the links form a staggered structure, 
doubling the size of the unit cell, and introducing a scat- 
tering between states of momenta k and k + t:. To deal 
with this we introduce the same two-component wave- 



functions as was used in Section VI 



Using this model we have investigated how the optimal 
stripe evolves as the second-neighbor hopping amplitude 
t' is changed from zero. For small values of t\ approx- 
imately, —0.3 < t' < 0.15, the structure of the stripe is 
largely unchanged. It is still described hy qx = and 
vanishing ^. The optimal number of holes per unit stripe 
length also remains {S ~ 0.46). All that happens is ba- 
sically that there are small redistributions of the holes 
within the stripe. 

Concerning the tendency to global phase-separation we 
have considered the behavior of Edv,{t')/Eps{t'), i.e. the 
ratio of the energy per hole of the domain wall and phase- 
separated states, respectively. We find that this ratio de- 
creases as we increase t' slightly from zero, and decreases 
when f becomes negative. This indicates that a negative 
t' favors the domain wall configurations compared to the 
global phase-separation but, at least for small t' , the do- 
main walls are still unstable against phase-separation. 

When t' becomes larger than 0.15, the holes diffuse 
into the antiferromagnet and widens the domain wall. 
The optimal number of holes per unit stripe length in- 
creases. The spin-twist ^ becomes non-zero making the 
antiphase boundary wider. Since we are working with a 
finite width w of the domain wall, we get problems when 
the holes start diffusing away from the center of the do- 
main wall. For our model to be valid, we must require 
that the domain wall is Neel-ordered and undoped at its 
vertical edges. This is to avoid surface energies when 
gluing together the domain wall with a Neel region. 

In the case relevant for the hole-doped cuprate planes, 
i.e. t' < 0, we find that decreasing the value of t' be- 
low —0.3 keeps the structure of the stripe rather intact 
in the sense that ^ remains vanishingly small and that 
the holes are localized close to the antiphase boundary. 
The optimal number of holes per unit stripe length does 
however change, it is reduced as t' is decreased, e.g. at 
t' — —0.5 we find that the energy per hole is minimized 
by (5 ~ 0.38. However, if we further decrease t' the opti- 
mal doping 6 will rise again as we reach t' ~ —0.7. For 
such large negative values of t' , the holes will spread into 
the antiferromagnet just as we found in the positive f 
case. Our numerical calculations indicate that ^ remains 
small, i.e. the spin-twist still occurs over a very short 
distance. 



VIII. CONCLUSION 



We have considered an effective version of the two di- 
mensional t — J-model where the electrons are considered 
as completely spin-polarized. The effect of such a spin 
texture is to generate a fictitious topological flux through 
the lattice. In this piaper we have extended the discussion 
of a previous paperld concerning the properties of these 
fluxes. 

Keeping in mind the result of Hasegawa et aL§, where 
it is shown that the energy of a free electron gas on a 
square lattice is minimized when there is one flux quanta 
per particle, the possibility of the system prefering a flux- 
generating spin-configuration does not seem remote. 

To check the above theory we have performed Hartree- 
Fock mean-field theory calculations for the system. In the 
Hartree approximation it seems like the system prefers to 
form a fiux-phase for certain choices of doping and cou- 
pling constants. However when the exchange-terms are 
included this effect seems to vanish and the coplanar spi- 
ral phase is energetically more favorable than the flux 
phase. Introducing a nearest-neighbor Coulomb repul- 
sion it is possible to make the flux-phase energetically 
most favorable also with the exchange terms present. 
However, the calculations indicate that for a wide range 
of dopings, these uniform phases are unstable against 
phase-separation into an undoped antiferromagnet and 
a highly doped coplanar spiral phase. Thus we have to 
conclude that from a thermodynamic point of view, we 
do not expect to find a fiux-phase in the phase-diagram 
of the model we have considered. 

The main part of the paper has been concerned with 
the generalization of this construction to describe stripes, 
directed along either the (10)- or (ll)-direction of the 
lattice. These stripes are appealing as they provide a 
smooth realization of an antiphase domain wall, contin- 
uously merging two Neel-ordered regions with opposite 
signs of the order-parameter. The holes naturally reside 
within this domain-wall which generates a fictitious flux 
which further can reduce the energy of the holes in do- 
main wall. 

Using J = 0.40, we find that the optimal antiphase 
domain wall is site-centered and runs along the (10)- 
direction of the lattice. The structure of this domain wall 
is such that there is a sharp spin- twist, basically only af- 
fecting a single row of spins, which aligns the spins ferro- 
magnetically in a one-dimensional channel. The doping 
of the domain wall is approximately (5 = 1/2 holes per 
unit length of the stripe. All holes are tightly bound 
to the domain wall, spreading over three rows of lattice 
sites. Due to the narrow spin-structure of this stripe, 
we found that the system does not exploit the 7r-fiux 
generated through each plaquette by this domain wall 
configuration. 
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We also find that stripes directed along the (11)- 
direction are energetically very close the optimal (10)- 
stripe, and if we decrease J they will become the most 
advantageous domain walls at J ~ 0.36, while increasing 
J favors the (lO)-stripe. Looking at the structure of the 
optimal (ll)-stripes we found that they do not generate 
any fictitious flux. 

Comparing the energetics of the domain walls with 
global phase-separation, we find that within our approxi- 
mation the global phase-separation is favorable. Further- 
more, we have incorporated second-neighbor hopping in 
the case of (lO)-stripes. For small values of this hopping, 
the structure of the optimal domain wall remains, while 
at larger values the holes starts spreading out, widening 
the domain wall. 

There are a number of questions that are left unan- 
swered at this point, and which we believe are interest- 
ing to further investigate. Concerning the next-nearest- 
neighbor hopping it would be interesting to investigate 
more through how it affects the tendency towards global 
phase-separation, and if it can stabilize the stripes. Fur- 
thermore, it would be interesting to further examine the 
effect of this hopping on the structure of the stripes, and 
also investigate its effects on the (ll)-stripes. 
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